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An improved method for extracting infrared exponents from functional equations is presented. The 
generalizations introduced allow for an analysis of quite complicated systems such as Yang-Mills 
theory in the maximally Abelian gauge. Assuming the absence of cancellations in the appropriately 
renormalized integrals the only consistent scaling solution yields an infrared enhanced diagonal 
gluon propagator in support of the Abelian dominance hypothesis. This is explicitly shown for 
SU{2) and subsequently verified for SU{N), where additional interactions exist. We also derive the 
most infrared divergent scaling solution possible for vertex functions in terms of the propagators' 
infrared exponents. We provide general conditions for the existence of a scaling solution for a given 
system and comment on the cases of linear covariant gauges and ghost anti-ghost symmetric gauges. 
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INTRODUCTION 



Functional approaches like Dyson-Schwinger equations (DSEs), see e. g. [Illl], or renormalization group equations 
''pi (RGEs) |4j have become standard tools to investigate the infrared (IR) behavior of Landau gauge Yang-Mills theory 
' [SHZ] during the last few years. The use of stochastic quantization [MIU] leads to equivalent results, especially the 
same IR behavior is found for a power law type solution in the deep IR. The IR enhanced quantities that lead to 
confinement of gluons and quarks have been identified as the ghost propagator jHl 13 HH 112] and in turn the quark- 
'— ' gluon vertex |131 114|. There are several reasons why Landau gauge was the first choice for such methods: First of 
all the ghost-gluon vertex does not acquire an IR enhanced dressing as suggested by Taylor [T5] and confirmed later 
on by DSE calculations (TBHIS] and lattice simulations [20 22 . This provided for a long time the starting point of 
the IR analysis. Second, Landau gauge is the simplest gauge in terms of DSEs, since it has only two fields and three 
[~ — interactions, which are the lowest numbers possible for a local gauge condition. Other gauges such as Coulomb gauge, 
00 the maximally Abelian gauge (MAG) or ghost anti-ghost symmetric gauges involve more degrees of freedom and more 
interactions rendering the functional equations more intricate. 

Numerical solutions for the coupled system of integral equations of the ghost and gluon propagators were first 
obtained in ref. [6 and steadily developed by including quarks |23] or improving the truncations to get better 
On transversality properties of the gluon propagator [24 . The quark-gluon vertex and the quark propagator have been 
investigated as a coupled system [l^, , whereas the other three-point functions [16. ,25) and the four-gluon vertex [2^ 
^ have been treated in a semi-perturbative approach. Again this is only possible because of the ghost-gluon vertex, which 
has a constant dressing function in the IR, and the dominant role of the ghosts. Having first results a posteriori enabled 
the investigation of the IR regime in a more general way. For not primitively divergent vertex functions a qualitative 
^ solution can be derived in d dimensions in terms of power laws for the dressing functions |27fE5] . The combination 
of DSEs and RGEs provides the means to identify a unique IR scaling solution |18) I19|. as does the assumption of a 
stable skeleton expansion [17]. Indeed, as we will show below, the latter is a consequence of constraints derived from 
RGEs. In the light of the success of a scaling analysis of the low momentum behavior of Yang-Mills theory in Landau 
gauge, it seems promising to pursue this method also in other gauges like the MAG. 

In general DSEs and RGEs allow two different types of solutions in Landau gauge |17 | f24 j . called the scaling solution, 
in which the dressing functions are related via scaling relations [SJ |S1 IH] , and the decoupling solution [171 |W '32] . 
The latter features a constant gluon propagator at zero momentum (sometimes interpreted as a gluon screening mass) 
and a tree-level like ghost propagator. Both solutions can be obtained within the framework of DSEs and RGEs, 
depending on the choice of boundary conditions [24_ , and it has been shown in ref. [175 that only the scaling solution 
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leads to IR enhanced vertex functions. However, in order to have static quark confinement general bounds on the 
IR exponents of the ghost and gluon propagators must be satisfied Both solutions fulfill this criterion. Another 
issue of the decoupHng solution is the breaking of global gauge and BRST symmetries as discussed in ref. [21] • Lattice 
simulations and the refined Gribov-Zwanziger framework |32l HO] seem to favor the decoupling solution, but 

for the former there are still unresolved issues concerning discretization jjl] and gauge fixing j42j.^ 

Whereas the Landau gauge has been used intensively for the last ten years in the context of non-perturbative 
investigations, the MAG was less common. The reason might be its intrinsic complexity compared to the Landau 
gauge. However, it offers an appealing picture of confinement, in which chromoelectric charges are confined by the 
condensation of chromomagnetic monopoles |44l Thereby the chromoelectric flux is squeezed into flux tubes by 
the Meissner effect. At the ends of these tubes the chromoelectric charges are located. This conflnement scenario 
works in analogy to the conventional superconductor, but with the role of electric and magnetic charges exchanged. 
Thus it is called a dual superconductor. In compact U{1) lattice gauge theory it is known that at low values of /3, 
i. e. for large coupling, magnetic monopoles indeed do condense and confine electric charges |451 ITZ]. One can try 
and generalize this to non-Abelian gauge theories by a specific choice of gauge, so that the theory can be expressed 
as an Abelian theory with magnetic monopoles [IHI- A special choice of this class of gauges is the MAG, where the 
off-diagonal components of the gauge-field are gauge-fixed such that their norm is minimized with respect to gauge 
transformations [15]. The chromomagnetic monopoles then can confine all chromoelectric charges if they condense. 
The appearance of these topological objects was first investigated by lattice calculations in refs. |481 149]. 

The existence of monopoles and their role for confinement was the topic of a recent lattice calculation, where no 
gauge was fixed [50]. The string tension calculated from the full gauge configurations, the Abelian part and the 
monopole part showed an amazing agreement. This hints at the importance of the Abelian part for confinement. 
Since monopoles are gauge dependent objects, it is not completely clear what role they play in different gauges, as 
for example they don't exist in Landau gauge. But it seems that the MAG is very well suited to identify the physical 
contributions of monopoles. Since it is easier to extract monopoles in the MAG than in gauge unfixed calculations 
like in ref. [50], it can provide a useful means of extracting the physical part of monopoles on the lattice (50]. On 
the other hand gauge fixing is necessary for all functional approaches. Therefore the MAG may be the best way to 
investigate monopoles with these methods. However, to establish a connection to monopoles is beyond the scope of 
this article. 

The dual superconductor picture of confinement leads to the hypothesis of Abelian dominance [5T: Since the 
classical configurations of the monopoles live in the Cartan subalgebra, i. e. the maximal Abelian subalgebra, of the 
gauge group, the relevant degrees of freedom should be Abelian. This can be interpreted as a mass for the off-diagonal 
gluons [52| . so that they decouple below this momentum scale and do not influence the dynamics in the IR. Indeed 
lattice results |52H54| indicate that the off-diagonal gluon propagator is suppressed at low momenta compared to the 
diagonal propagator. An analytic study did also find a mass for the off-diagonal gluon [55]. This result holds, if one 
takes into account the existence of the Gribov horizon (Gribov-Zwanziger framework) jSBl 157]. The diagonal gluon 
propagator then shows a Gribov-type behavior, i. e. it vanishes at zero momentum. However, by introducing additional 
terms into the Lagrangian that correspond to non-vanishing condensates (refined Gribov-Zwanziger framework) all 
three propagators of the MAG were found to possess a massive IR behavior [58^ . 

In this article we mention another possible manifestation of Abelian dominance: The diagonal gluon propagator is 
enhanced in the IR similar to the ghost propagator in Landau gauge. Therefore at low momenta it constitutes the 
dominant degree of freedom, which determines the behavior of the other two propagators and all vertices. However, 
since our analysis can only be the first step towards a full solution of the propagators in the MAG, we only mention 
this possibility and leave further discussions for the future, when more detailed results are available. 

There have also been investigations of the MAG on the lattice |52fl51] . The data for the lowest reached momenta 
showed that all three propagators are IR finite |54| . This agrees with results from the refined Gribov-Zwanziger 
framework [S^, where the addition of condensates also leads to the same IR properties. However, there has not been a 
successful attempt to use dynamical functional approaches to assess the low momentum behavior yet.^ The question 
is if there is a similar discrepancy between different methods as in the Landau gauge and if yes, what qualitative 
properties a scaling solution in the MAG might have. To fully answer this requires a complete numerical solution of 
the DSEs (or other functional equations) at all scales. For such an endeavor knowledge about the asymptotic behavior 



^ Recently an interesting suggestion was made how to implement different choices of the boundary conditions similar to functional 

equations also on the lattice 1431 . By selecting the Gribov copies with respect to the value of the ghost dressing function, different 

behaviors for the ghost propagator were found. 
^ An earlier attempt to investigate the MAG with DSEs was performed in ref. 1591 , where all two-loop contributions were neglected 

and the issue of bare two-point functions in the DSEs was not addressed, see sec. |III| In this article we will show that under these 
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of Green functions proves useful |51|Bn]- For high momenta we can rely on perturbative calculations, which have been 
performed up to three loops in the MAG [61j . Since perturbation theory is an expansion in the coupling constant it is 
clear how to classify contributions by the number of loops and the UV leading part consists of one-loop diagrams. For 
asymptotically low momenta, however, no information is available. In this article we will attempt such an analysis, 
where one interesting result will be that the structure of IR dominant diagrams is quite different from all other known 
examples: Instead of one-loop diagrams (as in the Landau gauge) two-loop diagrams are IR leading. This constitutes 
the first known case of this type and prohibits the use of known approximations for a numerical solution. 

The investigations of the MAG mentioned above were all done for the gauge group SU{2). It is an interesting 
feature of this gauge that the number of interaction terms in the Lagrangian is different for SU{2) and higher SU (N). 
Thus here the question of the gauge group dependence is of special importance, in contrast to Landau gauge, where 
both functional approaches and lattice simulations |62fl51] agree that the qualitative behavior is the same for SU{2) 
and SU{3). Since the DSEs in the MAG are even more complex for SU{3) than for SU{2) the need for a more 
effective method to investigate the IR regime is obvious. The scheme presented below provides this and shows that 
the additional interactions of SU (3) neither alter the solution derived in SU (2) nor allow for an additional one. 

The article is structured as follows: In sec. |TT] we illustrate the basic idea of the method of power counting to 
determine qualitative IR scaling solutions. It serves as an entry point for the reader new to this kind of analysis and 



may be skipped by the reader familiar with this topic. Sec. Ill explains the method we will employ to investigate 
the MAG. It is not specific to the MAG and thus the results may be used for other gauges as well. However, we will 
repeatedly exemplify some steps in the Landau gauge to provide a connection from the abstract method to familiar 
facts and relations. The MAG and other gauges are investigated in sec. |IV| and sec. |V) respectively. Our conclusions 



are presented in sec. VI Some detailed calculations are deferred to three appendices. 



II. THE INFRARED POWER COUNTING TECHNIQUE 

In our study of the IR behavior of Yang-Mills theory the quantities of interest are the full non-perturbative Green 
functions. Usually one starts the analysis by considering two-point functions, i. e. propagators. They are parameter- 
ized by 

^viP)=P^,^^ (1) 

where Pij is the part containing color and Lorentz structure and Z is some dressing function. For low momenta, i. e. 
— >■ 0, the latter obeys due to renormalization group arguments a power law, 

Z'"ip')^A-{p^r, (2) 

with some constant A and the quantity a called infrared exponent (IRE). Clearly if a < the propagator is enhanced 
in the IR compared to the tree- level propagator and it vanishes for a > 1/2. We assumed here for simplicity that 
the propagator can be described with a single dressing function, but a generalization to several dressing functions is 
possible. Although vertices are more complicated and feature in general many dressing functions corresponding to 
different tensors, all of them show the same momentum dependence when all momenta vanish. This case is in the 
following denoted as the uniform limit and will be studied in this work. 

Although DSEs can be used to determine dressing functions for the complete momentum region, we are only 
interested here in the behavior at low external momenta of the integrals. Since they contain terms like l/{p — fc)^, 
where one momentum is the loop momentum and the other the external one, main contributions to the integrals arise 
from loop momenta that have the same order of magnitude as the external ones. When all momenta are small, we 
can replace the dressing functions by their respective IR expressions. Due to dimensional reasons all momenta of the 
integral transform to external momenta upon integration. For two-point functions this can directly be seen from the 
available analytical expression for the one-loop integral ^ 

Similar expressions can be obtained for the two-loop integrals of two-point functions by successive integration of the 
integrals. In case of higher vertex functions analytic results are only known for one-loop three-point functions |651I66| . 



From UV and IR convergence one can see that the exponents ui and U2 have to fulfill d/2 + i/i + 1/2 ^ 0, d/2 + i/i > and d/2 + 1^2 > 0. 
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An explicit calculation of three-point functions in Landau gauge has been presented in refs. flS'. 'WP. A very successful 
method for the evaluation of these kind of integrals is the negative dimensions integration method (NDIM) |68-n71j . 
In this work we only treat the case of uniform IR divergences, i. e. all external momenta tend to zero simultaneously. 
The case of only a subset of the external momenta going to zero is discussed in refs. |171 [^ . For an explanation of 
how the power counting works for RGEs see ref. |18j . 

Knowing that all integrals transform into powers of the external momenta one can work directly at the level of the 
IREs (but has then to ignore possible cancellations of integrals). We want to make this point clear with a simple 
example. Consider the gluon propagator DSE in Landau gauge, which is given diagrammatically in fig.[lj Evaluation 






FIG. 1: The gluon propagator DSE in Landau gauge: The bare propagator, the ghost-loop, the gluon-loop, the tadpole, the 
squint and the sunset diagrams. Small dots represent bare vertices and large dots IPI vertices. 



of integrals yields 



2^1+2(53g+4(5g; j^squint 



(4) 



with the IREs 5gi, Sgh, S^g, Sgg and S^g of the gluon and ghost propagators, the three-gluon, the ghost-gluon and 
the four-gluon vertices. The Ls denote constant terms for the gluon-loop, the ghost-loop, the sunset and the squint 
diagrams.* They depend on the IREs (see eq. ^ above) and contain also the constant coefficients from the power 
laws. When p is small not necessarily all terms on the right-hand side of the equation scale equally, but at least one 
must have the same exponent as the left-side, what can be phrased as 

{p^)-'^' = max ((p2)0, (p2)53,+25,, ^ (p2)i,,+25,, ^ (^2)54^+35,, (^2) 2^3, +45,,^ _ (5) 

Canonical dimensions do not occur here any longer, since they cancel in four dimensions. Eq. (|5| is equivalent to 

-Sgl = min(0, Ssg + 2Sgl,Sgg + 2(5gft,, S^g + '3Sgl,2S3g + 4(5^). (6) 

From this equation we can extract single inequalities, e. g. 

0<S3g+iSgl, (7) 

which can be used to show that 6gi > 0: Since the three-gluon vertex is primitively divergent its DSE features the 
bare vertex, so that we have < 0, i. e. the vertex' IRE is non-positive. Taking this into account we get 



< -Sag < S, 



Ogl 



6gl > 0. 



(8) 



In sec. Ill we will show how to single out the relevant inequalities, so that one does not have to work with inconvenient 
minimum functions. 

First IR analyses of DSEs and RGEs in Landau gauge relied on investigating individual diagrams. Whereas at first 
only the IR behavior of the propagators and the ghost-gluon vertex was determined jS] [51 [TT] , this could be extended 
later on to the whole tower of Green functions by use of a so-called skeleton expansion |27] . Thereby higher n-point 
functions are expanded in loops consisting of dressed primitively divergent vertex functions. Another approach was 



* The tadpole was not taken into account as it does not depend on the external momentum and can thus be absorbed in the renormalization. 
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taken in [T51 [TO] , where the combination of the two equivalent systems of DSEs and RGEs allowed to show the 
uniqueness of the solution for the IREs. In |T7| it was shown that also the assumption of a non-explicitly divergent 
skeleton expansion allows a proof of uniqueness. Non-explicitly divergent means here that insertions, necessary to get 
higher orders of the expansion, do not have negative IREs. Since, as shown below, the constraints obtained from this 
assumption follow also from the RGEs without such an assumption, this result is not surprising. Furthermore in [17] 
the existence of additional IREs for certain kinematic limits, e. g. only one momentum of a three-point function going 
to zero, was also established. A numerical confirmation of kinematic IR singularities for the three-point functions of 
Landau gauge Yang-Mills theory has been presented in ref . ^25j , where it was also shown that they only appear in the 
longitudinal parts of these vertices. Also the quark-gluon vertex features kinematic singularities |14l 17^ , which turn 
out to be as strong as the uniform ones. This is in contrast to the pure Yang-Mills case, where they are significantly 
weaker [1^1 E5]. 

The first step towards the improved technique used in this article can be found in [28 , where a formula for the 
IREs of all diagrams in a DSE was presented. Working in Landau gauge the known exponents of propagators and 
the ghost-gluon vertex served as input. However, that result still was based on the use of a skeleton expansion, i. e. 
only three- and four-gluon vertices and no other vertex functions were taken into account. A generalization of this 
formula to arbitrary interactions and without assumptions about any of them is used in the next section to derive a 
lower bound for vertex IREs. 



III. MAXIMALLY INFRARED DIVERGENT SOLUTION AND SCALING RELATIONS 

In this section we give an overview of the method we are going to employ in sec. |IV| to derive an IR solution 
for the MAG. We specialize on pure scaling solutions and do not consider dynamical mass generation or kinematic 
divergences. The analysis is similar to the one presented in the recent ref. jTH] for Landau gauge. We will keep the 
discussion general without fixing the number or type of fields yet, since doing so does not provide any simplification; 
on the contrary, writing some equations explicitly for the MAG leads to unnecessarily long expressions. The crucial 
results of this section might prove useful also in the investigation of other actions, as we identify the most restrictive 
inequalities one can derive from the system of functional equations in the pure scaling case. Up to now the systems 
of equations were studied diagram by diagram, whereas here we provide a more general approach that enables one to 
write down the decisive inequalities. These contain the only information as extracted from functional equations that 
restrict the IR behavior. All other inequalities, of which there are infinitely many, arising from individual diagrams 
are superfiuous in any investigation and can be disregarded. As we do not employ any truncation, we analyze the 
full system of DSEs and RGEs. In principle it would not even be necessary to write down the equations explicitly, 
since the necessary information is encoded in the interaction types of the action. Thus with the results of this section 
one can then get the IR scaling solution rather easily only by considering the type of interactions in the Lagrangian 
without explicit power counting of individual graphs. 

Although the results are simple, the derivation is somewhat technical. Below we try to explain the idea and defer 
some details to appendices. To allow a more intuitive understanding of some relations we give examples in the Landau 
gauge, since many readers might be familiar with it. We start with a general formula for the IRE of an arbitrary 
diagram derived purely by combinatorics. It will lead to the notion of the maximally IR divergent solution. Then we 
analyze the propagator equations and show that at least one primitively divergent vertex is not IR enhanced. 

For the existence of a scaling solution we have to make some assumptions: 

• An IR enhanced propagator, i.e. S < 0, is only possible if the term arising from the bare propagator is canceled 
in the DSE. This is simply a consequence of the general form of propagator DSEs, where the inverse dressed 
propagator is on the left-hand side and the inverse bare propagator on the right hand side, yielding the equation 

-6 = min{0, . . .) ^ 6>0. (9) 

Thus the bare propagator has to be canceled somehow to allow for IR enhancement. This is realized in Landau 
gauge for the ghost propagator by the so-called horizon condition [73 , which is implemented into the DSE of 
the ghost propagator by an appropriate renormalization condition. We work here under the assumption that 
similar mechanisms exist in other gauges, so that bare propagators can vanish effectively from the DSE at zero 
momentum. We comment on the possibility S — below at the end of subsec. |IIIB| 

• Propagators with several dressing functions lead to several DSEs that are obtained by projecting out the cor- 
responding tensor parts from the original DSE. Thus the individual parts of the propagators can be taken into 
account in a power counting analysis by different IREs. The vertices are split accordingly. However, different 
IREs are only expected if the structure of the inequalities changes. This could happen for example if the bare 
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vertex vanishes in the projection of a vertex equation. On the other hand, an IRE can also be determined by 
constraints like symmetries. This happens for example in linear covariant gauges j]p, where the gauge symmetry 
determines the longitudinal part of the gluon propagator (see sec. |V|. This leads to additional information that 
can be used in the analysis. 

• Cancellations always are a possible threat for a scaling analysis if terms are taken into account that actually 
disappear. We disregard this possibility here and assume that none of the important contributions behaves 
in this way. Indeed the analysis is only rendered invalid if the IR leading terms vanish, while a vanishing of 
subleading terms is not relevant. 



A. A formula for the infrared exponent of an arbitrary diagram 



Before we derive any formulae we explain our notation. We do not specify any fields explicitly in our derivations 
for two reasons: First, it is not easy to write down the derivations in the MAG, because this gauge has so many 
interactions between its three fields that the equations become quite large. In our generic notation the derivation can 
be done in a compact way. Second, proceeding in a general way we can learn more about the basic structure of IR 
scaling solutions and the results can be used for various actions. So we denote the fields by the set {(j)s}, s = 1, . . . , S, 



with S the number of different fields in the action. In the MAG this is {A,B,c^c} (see sec. IV for details) and in 
Landau gauge {A,c,c}. Furthermore we need the following quantities: 

• 6i: The IRE of the propagator of the field 0^. 

• Si...k- The IRE of a vertex with legs of the fields (j>i, . . ., 

• n-i,,,k'- The number of vertices of the type 0i . . . (j)k- Superscripts d or 6 denote only dressed or bare vertices, 
respectively, whereas no superscript refers to both. 

• nf. The number of internal propagators of the field 1/)^. 

• nii: The number of external legs of the field (f>i. 

• Ci,,,j: The canonical dimension of the vertex indicated by the subscript. 

• kf'"'': The number of legs of the field the vertex (pj . . . (f>i has. 

One of the main pieces in the forthcoming analysis will be a formula for the IRE Sy of an arbitrary diagram v with 
canonical dimension Cy. With the notation introduced above we can write down this general formula for an arbitrary 
Z-loop diagram in d dimensions: 

Sy=d+Y,MS^-^)+ J2 <...^.('^'l■■.V+C^l...^„)+ (10) 

i vertices,r'>3 

+ -c„. (11) 

vertices,r>3 

The double subscripts appearing in the last two sums indicate all possible combinations of r fields. E. g. for the 
Landau gauge the term corresponding to r = 3 is 



n 



d 

AAA 



( 5 AAA + \]+ "Icc \5a-cc + \]- (12) 



Using topological relations this expression is reformulated such that the dependence on internal propagators is ex- 
changed for one on the number of external legs. The details of this calculation are given in Appendix The final 
result in four dimensions is 



vertices, r'>3 \ i / vertices, r'>3 



(13) 



The sums extend over all vertices with r legs in the diagram. Eq. ( 13 1 is a purely combinatoric result without any 
assumptions whatsoever employed. At this point one may be tempted to employ a skeleton expansion, i. e. restrict 
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the sums to primitively divergent vertices. This is indeed possible since the derivation of eq. (13 1 was valid for a 
generic diagram, i. e. also those one gets in a skeleton expansion. Nevertheless an analysis without skeleton expansion 
can be done and we do not have to rely on the assumption of a stable skeleton expansion here. Note that for RGEs 
— and for every diagram in a DSE there is only one n!' — 1, while the other are 0. For a consistent solution of 



DSEs and RGEs eq. ( 13 1 hints already at the condition that at least one coefficient of an is zero, leading in turn 
to a non-scaling vertex. This consistency condition was first stated in ref. [18_. Below we will show in subsec. IIIB 
in detail how it emerges for general structures of interactions. 

Eq. (13) leads to the notion of a maximally IR divergent solution. Therefor we need information about the 
coefficients of n'^, which we can get from RGEs, and n'' in the sums. We explain the case of three-point functions and 
shift the derivation for higher vertex functions to Appendix [B] Consider the RGE for a generic three-point function 
with the fields A, B and C, where it is not excluded that some of them are the same. For instance A may be a gluon 
field and B and C the Faddeev-Popov ghost and anti-ghost fields of Landau gauge. In fig. |2] we only show the one 
diagram on the right-hand side we are interested in. At the level of IREs we have on the left-hand side only the IRE 
of the three-point function; on the right-hand side appears three times the same three-point function and one IRE for 
each propagator. Since the diagram on the right-hand side cannot be more IR divergent than the left-hand side, we 
can write down an inequality for the IREs: 



5 ABC < 36abc + Sa + 5b + Sc 



Sabc + 1:{^^ + ^b+Sc)>0- 



(14) 



Note that in four dimensions for uniform scaling the canonical dimensions always cancel and need not be considered 
explicitly, but the argument is also valid in d dimensions; for details see Appendix [B| Eq. (14 1 implies that the 
coefficient of n'^ABC ' -'^'^ ' non-negative for arbitrary fields A, B and C . If we can find other similar constraints 



for higher vertex functions, we know that the sum over dressed vertices can only increase but never lower the IRE of 
a diagram. Indeed such inequalities exist and their derivation can be found in Appendix |B] The general form is 



'^5,. > 0. 



(15) 




FIG. 2: One specific diagram in the RGE of a generic three- 
point function. Internal lines represent dressed propagators, 
black blobs dressed vertices. The gray blob is a regulator. 

Interestingly it turns out that for primitively divergent vertex functions the same restrictive inequalities arise from 
the assumption of a non-divergent skeleton expansion or, phrased differently, the skeleton expansion cannot lead to 
more divergent terms. Such an expansion expresses non-primitively divergent vertex functions by a loop expansion 
containing only dressed quantities. Higher orders can be obtained by inserting additional loops using basic insertions, 
which consist only of primitively divergent vertex functions, into a diagram [27J. By non-divergence we mean that 
the IRE of such an insertion is greater or equal to zero so that the new diagram is not more IR divergent than the 
original one. Using the skeleton expansion one gets a series of infinitely many diagrams and it is a priori not clear, 
whether the sum can change the IR scaling. The general result from above shows that this is actually not the case 
and previous results based on the skeleton expansion are indeed valid. Note that at least in Landau gauge it was 
found for two, three and four dimensions that the insertions do not change the IREs, i. e. the inequalities above are 
saturated |27ll28]. 

The inequalities above also restrict the coefficients of in the second sum of eq. ( 13 1. We know that for primitively 



divergent vertex functions the bare vertex appears on the right-hand side of its DSE. Therefore these n-point functions 
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have an IRE lesser than or equal to zero, 



^prim. div. vertices 



<o, 



(16) 



and not only the coefhcients of the n'^ are non-negative, but also those of the n'', since they only differ by a non-positive 
number. The fact that eq. ( 16 I is only valid for primitively divergent vertex functions causes no problem, since 



only for primitively divergent vertex functions. We demonstrate explicitly how this argument works for a three-point 
function: 



1 



{Sa + Sb + ^c) <Sabc <0 



5a + 5b + 5c> 0. 



(17) 



For a general primitively divergent interaction 



this can be written as 



(18) 



This inequality can also be derived with the same argument from DSEs, since DSEs yield similar inequalities as those 
in eq. (15 1, but with a different numerical factor in front of the vertex IRE. 



Thus the second and third terms in eq. (13 1 are both non-negative and the maximally IR divergent solution, i. e. 



that with the lowest IRE possible, for n-point functions with rrii legs of the field 0^ in terms of the propagator IREs 
5i is 




(19) 



It is realized when the inequalities connecting the IREs of vertices and propagators, eq. ( |15[ ), are saturated and thus all 
coefficients of the n'^ vanish. Then all diagrams in an RGE scale equally (since n'' — 0). Diagrams in DSEs, however, 
contain a bare vertex, which is reflected by the appearance of n''. Its coefficient vanishes only for the IR leading 
diagrams, e.g. those with a bare ghost-gluon vertex in Landau gauge, while the other diagrams have larger IREs. In 
Landau gauge these are the diagrams with bare three-gluon and four-gluon vertices, which have an IRE that is 3k 
and 4k higher than for the leading diagrams, respectively. These numbers stem directly from the coefficients of the 
n**. Thus for the maximally IR divergent solution the bare vertex divides the diagrams into classes of IR importance. 



The important formulae of this subsection are eqs. (13 1, (151 and (18 1. The latter two comprise all constraints on 



the IREs one can get in a uniform scaling analysis. They can directly be derived from the vertices of the theory and 
there is no longer the need to write down many inequalities from all sorts of different diagrams. To demonstrate this 
we again employ the Landau gauge as an example, where one gets up to four-point functions the following inequalities: 



34g 



5gl >0, 

25gi >0, 



2<5<,h + 6gi >0 
Ssg + ^-5,,, >0, 



Sgg+Sgh + ^Sgl >0, 



^2gh2g 



Sal >0, 



(20) 



with 52gh2g the IRE of the ghost-gluon scattering kernel. Any attempt to obtain more restrictive inequalities fails. 
Thus on the level of the IREs one can reduce quite large systems of DSEs to a much smaller set of inequalities without 
making any assumptions or missing any information relevant for the uniform IR scaling. 



B. Analysis of the propagator equations 

In this subsection we will have a closer look at the propagator equations. For a given system of DSEs one could 
of course perform a one-by-one analysis of each diagram. This is a simple task in Landau gauge, where the ghost 
and gluon propagator equations have one and four diagrams, respectively, neglecting the bare two-point functions 
and tadpole diagrams. The situation is more involved for the MAG, where in SU{2) we have six, seven and seven 
diagrams for the diagonal gluon, the off-diagonal gluon and the ghost propagator, respectively. These numbers increase 
for S'C/(3) to 13, 17 and 14. Therefore such an analysis easily becomes extensive. We describe below how the analysis 



9 



can be performed efficiently for all diagrams at once, again by using abstract notation. As a result we can show that 
for a consistent scaling solution at least one of the inequalities found above has to be saturated and thus at least one 
vertex cannot be IR enhanced. 

In the case that we only have three- and four-point interactions there are only four different types of diagrams 
besides the bare two-point function: the tadpole, the one-loop, the sunset and the squint diagrams, see fig. [3] All of 
these possibilities appear in the DSEs of the MAG in at least one realization, see fig. [4] However, the argument given 
below holds also for higher interactions. 



(a) 



ib) 






FIG. 3: Types of graphs in a two-point DSE: (a) tadpole, (b) one-loop, (c) sunset, (d) squint. 



The analysis is based on the fact that the IRE of the left-hand side has to be matched by at least one IRE on 
the right-hand side. Thus we do not have an inequality at this point but an exact equation, but for now we do not 
specify which diagram the IR leading one might be. The IRE of a diagram can be written as the sum over the IREs 
of all propagators and all dressed vertices contained, as was done in sec. |III A[ The equation can be rewritten using 
topological relations and transformed into an inequality using the lower bound on IREs given by eq. ( 19 1. The general 
form is: 



1 b ^ 



> 0. 



(21) 



The derivation of this inequality is presented in Appendix |C] This inequality has decisive differences to the previous 
ones: The sign is opposite to that in eq. (18 1 and it depends on the bare vertex via nj^ j^. While inequality (18 1 
is always valid, inequality (21) is only true for the leading diagram, which is not determined up to now. The factor 



divides the large number of possibly leading terms into classes determined by the bare vertex of each diagram 



as discussed above. For the leading class one can combine eqs. (18 1 and (21 1 to get 



0, 



(22) 



where one has to keep in mind that the bare vertex determines the appearing IREs. To clarify that this equation 
is only valid for the leading class we added a hat to n^. This sum is nothing else than adding up the IREs of the 
propagators corresponding to the legs of a vertex. From eq. ( 22 1 it also follows that the corresponding vertex has a 
zero IRE, since inequality ( 15 1 reduces to 



(23) 



Together with the non-positivity of the IREs of primitively divergent vertices, eq. (16 1, this yields 



Sii...i^ = 0. 



(24) 



All one has to do now is to test all possible equalities one can get for a consistent solution. The number of equalities 
is the same or less than that of the bare vertices of the theory. This is a manageable set of equations that describes 
the tower of infinitely many coupled DSEs and RGEs. 

If the signs of some propagator IREs are known, the number of possible solutions reduces considerably. One can 
also see that the identity eq. (22 1, called scaling relation, can only be fulfilled if a positive and a non-positive IRE 
are combined, i. e. having only positive propagator IREs leads to the only solution that all IREs are zero. This 
corresponds to the perturbative behavior of propagators and vertices in the ultraviolet and seems unlikely to be 
realized also for small momenta, since perturbation theory itself predicts its breakdown in this regime. 
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Let us illustrate the efficiency of the method by analyzing the Landau gauge. Here there exist only three different 
bare vertices, which yield two different equalities: 

Sgl^O, hgl+Sgh = 0. (25) 

The first equation leads to the trivial solution Sgi = 5gh = 0, while the second one yields the known IR scaling 
solution 5g + 25 gh = [5] with a finite IR dressing function for the ghost-gluon vertex. The IR leading diagrams in the 
propagator DSEs are those with a bare ghost-gluon vertex and the solution for higher vertex functions |19 | l?7 l l^i can 
directly be derived by adding to these diagrams further legs using ghost-gluon vertices. This example should make 
clear that with the method developed here a complete analysis, i. e. including higher vertex functions and two-loop 
diagrams, can be done in an easy and elegant way without relying on a skeleton expansion or assumptions on the 
ghost-gluon vertex. 

If there are several diagrams with the same bare vertex in the propagator DSEs it may be that more than one 
solution emerges. The reason is that our analysis assumed that one specific diagram is IR leading, but the end result 
only depends on the bare vertices. Thus one has to distinguish between the results obtained from different diagrams. 
Note that the obtained scaling relation is not influenced by this, but the behavior of higher vertex functions. Take as 
an example sunset and squint diagrams. Assuming the squint diagram to be leading, we get that four-point functions 
are not IR enhanced and thus the sunset is also IR leading. On the other hand starting with the sunset, we do not get 
information about three-point functions, which are the only dressed vertices in the squint. So it is not clear, whether 
the squint is IR leading or not in this case. This is exactly what happens in the MAG and consequences are discussed 
in subsec. IIVBI 

One of the main results of this section is given by eq. ( p2| . In order to get a scaling solution, we write down 
all realizations of this equation, which are determined by the bare vertices of the theory. Then we test each single 
equation, if it is consistent with the inequalities derived above. If we find such an equation, it is called a scaling 
relation for the IREs of the propagators. The vertex corresponding to this equation is called the leading vertex and 
does not get IR enhanced. The equation may have been derived from several diagrams in the propagator DSEs. If 
this is the case, details of the solution may depend on it. Indeed the MAG is an example where this happens and 
details can be found below in subsec. |IV B[ The IREs of higher vertex functions are derived by adding further legs to 
the leading diagram(s) in the propagator DSEs by insertions of the leading vertex. 

IV. INFRARED SCALING SOLUTION FOR THE MAG 

Before we employ the method presented in the previous section we will give some details on the MAG in the 
following subsection. The case of SU{2) is treated then with the method presented above, whereas the inclusion of 
the additional interactions for higher SU{N) is deferred to an own subsection. Finally we comment on the treatment 
of the longitudinal part of the propagator for the off-diagonal gluon. 

A. The maximally Abelian gauge 

The crucial point of the MAG is the separate treatment of different components of the gauge field. By splitting it 
up into 

Af, = Al^T'' = B'^T" + A^^T\ (26) 

where the are the — 1 hermitian generators of the SU{N) Cartan subalgebra and the T° the N'^ — N remaining 
generators, we divide it into the so-called diagonal and off-diagonal gluons A and B, respectively. This name is based 
on the possibility to write the generators of the Cartan subalgebra in diagonal form, 

/ 2 \'/' 

= L(j^i) ) X dzag(l_^,-j,0,...,0), j = l,...,N-l, (27) 

j times 

whereas the are off-diagonal matrices. Indices are used in accordance with recent literature, i. e. a,b, . . . for 
off-diagonal components and . . . for diagonal ones; see e. g. |58l r,s, . . . refer to both types. We choose 
Hermitian generators: 

[T'',T']=ifT* (28) 
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with the orthogonahty condition (T/ = 1/2 for SU{N)) 

tr{T''T''} ^TfS''''. (29) 



The underlying idea of the maximally AbeHan gauge fixing is to minimize the off-diagonal components |46[ |48| 149] 
and thereby to emphasize the role of the diagonal gluon by demanding that the functional 

R[B]= I dxB'/.B; (30) 

should be minimized with respect to gauge transformations. This yields the gauge fixing condition 

DfB^ = (31) 

for the off-diagonal gluons with the covariant derivative defined with respect to the diagonal gluon components: 

D';,' ■.= S-''d,+gr'^A;,. (32) 

For the residual gauge symmetry of the diagonal part we choose Landau gauge, i. e. 

d,A^^ = 0. (33) 

To keep the gauge renormalizable a quartic ghost interaction term has to be added to the Lagrangian |75| 175], All 
this can be done within the BRST formalism, so that the Lagrangian can be written in the convenient form |76| 

+ sc'^(7)fB^-z|6'') + s(-iA5r'''c'^cV- ^Agr^'^c'^cV), (34) 

with the BRST transformation s defined as 

s a; = -{D'^^'c'' - g r'^Blc^ - g r'^B^c'), s = ~{5^^d,c, - g r'B^c"), 

sc" = -^gr'"'c''c''+gr'"c''c\ sc^^ib", s6°=0, 

sd = -]^g r^c''c\ s(f = ih\ sb'^Q. (35) 

The parameters appearing in the Lagrangian are the gauge fixing parameter for the diagonal part ^ (for Landau gauge 
it is 0) and that for the off-diagonal part a. The parameter A controls the strength of the quartic ghost interaction. To 
show explicitly all interactions we give the full form for SU (N) , where we introduce a parameter C. that interpolates 
between the MAG — 1) and the Landau gauge (^ = 0): 



^ = ^id,B:)id,B^^) - -{d,Bl){d,Bl) -f -gr'^BlBl{{d,Bl) - (d^B^^)) + g B'^A^^Hd.B:) - (d^B^))- 



^-id,B:)id,B^^) - \{d,Bl){d,B;) -f i< 
+ r'^r^HBlAlB'^Ai - B^AlAl^Bi) + ^g^ r'^B^.B^M'"''' B'^AI + 2r''^Bf,Al - 2/''*A;i?^) + 

+ \{dMidM ~ \{d,Al){d.A^^) + \gr'B'^^Bl{{dM - {d,A^^)) -f \g' BlBlBlBt+ 

+ c°(,5'^^9^ + C9r'''Al){5''^d^ -f gf''A=^)c^ - g r^c^iiS'^d^ + Cg f'' A^,)B';,y- 
-gf^'c'^iiS^'d, + Qgr^'^ADBlc'^) - g^ Q ^'^^ B^^B^c'^ c^' - (1 - C)g r''c'^Bldy+ 

+ C'df,{df,c' - g P'^'Blc^) + ^g^ p'^'^C^C^c'^ + ^g^ r'"'f'"^'c''c''c'^c' + j-a&cjade-6-c^d^e_^ 

+ ^((^'''9, + Cgr''A^^)Bl)^ - ^g^r'T^'c^c^c^c^ - ^gr'\{6^'d^ + Cgr'Aj,)B'^^)c'>d- 

^gjabc^^^adQ^ -f Cg/''''M;)B^)c''c^ + ^g^r'f'^c'^c'c^c^ + Y^id^^^f . (36) 

The ghost field, being also in the adjoint representation, is split similar to the gluon field. However, the diagonal 
ghosts decouple for A = a, which we will use henceforth. This identity leads to the diagonal ghost equation |76| . 
which is an additional Ward identity respected by the action for A = a. 
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Three-Point 


Four-Point 


SU(2) 


ABB, Acc 


AABB, AAcc, BBBB, cccc, BBcc 


SU(N) 


ABB, Acc Bcc, BBB 


AABB, AAcc, BBBB, cccc, BBcc, ABcc, ABBB 



TABLE I: Interactions of the MAG fixed Lagrangian. There are 2 (4) three-point and 5 (7) four-point interactions for SU{2) 
and SU{N), respectively. A is the diagonal gluon field, B the off-diagonal one and c the ghost field. 



For SU{2) the structure constants assume the simpler form /"''^ — 0, /°^' — with e the antisymmetric tensor, 
because it is the only possible combination for one diagonal and two off-diagonal fields. For general SU{N) there 



do not exist non-zero structure constants with more than one diagonal index, as can be seen from eq. ( 28 1 . The 
Lagrangian contains all types of interactions possible for SU{N). Only in the case of = 2 those containing only 
off-diagonal indices vanish. The interactions are given in table |T] 

For the derivation of the DSEs we expand the action in the fields, which explicitly shows the types of interactions: 

S —-Sf^'^AiAj + -S^fP BaBb + S^lCaCb — -S^f^ AiBaBh — S^l'^AiCaCh — —S^i^^^ BaBijBc — S^^^BaCbC]:,+ 

1 qAABB a a u n ^ qAAcc a a t; ^ ^ qBBBB u U n U ^ qBBcc u U ^ ^ qcEcc- t; ^ ^ 
~ ^^ijab ^i^j^aJ^b ~ ^ ij'^i' ^i^j^a^b — J abed ^a^bJ^cJ^d ~ 2 abed J^a^bCcCd — abcd'^aCbCcCd — 

~ '^y^tahc^'^i^aBbBc — Sf^f^"^ AiBaCbCc- (37) 



The coefficients S'f'y'^^ correspond to the bare n-point functions. The DSEs were derived with the Mathematica 
package DoDSE [77] and are given in fig. |4] 

Before starting the IR analysis we have to clarify a detail in the DSE of the diagonal gluon propagator: If the bare 
propagator is not canceled and included in our analysis its IRE is non-negative. This leads to the trivial solution that 
all IRE are zero as explained below. Therefore we argue that for the existence of a scaling solution with IREs unequal 
to the propagator DSE of the diagonal gluon may be renormalized in a similar way as the ghost in Landau gauge, so 
that the bare propagator is absent in the scaling analysis. In Landau gauge this renormalization is connected with the 
so-called horizon condition, which takes into account the existence of the Gribov horizon [B HTUl [75H5T| . Recently it 
was shown in ref. [74^ that the Gribov region in the MAG is quite different from that in the Landau gauge: Whereas 
in the latter it is bounded in all directions of field space [101 152] . this is only the case for off-diagonal directions in the 
MAG [74] . The diagonal direction on the other hand is unbounded what can lead to a different low energy behavior as 
in Landau gauge. Another interesting fact is that purely Abelian configurations {B = 0) that are gauge transformed 
into Landau gauge lie at the Gribov horizon of Landau gauge |83| . But these configurations are exactly the ones 
important for the IR scaling solution in that gauge. 



B. Scaling solution for SU{2) 

With the tools developed above in sec. |III| the remaining analysis after having established the interactions of the 
MAG is done rather fast. Introducing the IRE of the diagonal gluon, the off-diagonal gluon and the ghost propagators 
as 5ai and dc, respectively, the occurring interactions give the following list of constraints from the inequalities 
(pl: 



-5a + 5b >0, 5a + 5b >0, 2,5b >0, 

\5a + 5c>0, 5a + S,>Q, 25,>Q, 

5b + >0. (38) 

Note that there is a _B c symmetry in SU{2), i. e. we could restrict ourselves to only one of these two fields and 
get the same results. We will not do so in the following. First we note that both 5b and 5c are non-negative^, which 
is a direct consequence of the self-interactions of these fields. This renders the inequality 5b + Sc > superfluous and 
allows to discard some of the inequalities, since they are less restrictive than others, e. g. ^5a + 5b > due to the 



^ The corresponding inequalities in d dimensions are 5s > (4 — c!)/4 and 5c > (4 — d)/4. 
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FIG. 4: Yang-Mills DSEs for the propagators in the maximally Abelian gauge. Continuous lines are diagonal gluons, dashed 
lines off-diagonal gluons and dotted lines ghosts. All internal propagators are dressed, the thick dots are IPX vertices and the 
small ones bare vertices. The unusual conventions for the lines are due to the use of DoDSE |77| . which cannot draw wavy 
lines. For the structure functions we did only take into account the antisymmetric /-function. Two diagrams (squints with 
a bare BBcc vertex) that are possible on the combinatoric level do not appear here, because they are zero due to the color 
algebra. 

existence of the AABB interaction that leads to Sa + ^b > 0. The remaining system of inequalities for the propagator 
IREs is 

5b >0, 6c >0, 

5a + 5b >0, 5a + 5c >0. (39) 

These inequalities always have to be fulfilled, but at least one of them has to be saturated to have a scaling solution. 
We examine every one of the four possibilities. 

The first \s 5b = 0. This renders 5a also non-negative by virtue of the third inequality. But with only non-negative 
propagator IREs left, there cannot be a non-trivial scaling solution for the remaining propagators, as this needs 
positive as well as negative IREs. It is a special feature of the MAG that every field type interacts with every other 
field type, so that whenever one IRE of a propagator is zero, all others become non-negative and only the trivial 
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solution with Sa = Sb = Sc = can be realized. Note that there are examples where one field has an IRE of but 
the IREs of the other fields still obey a scaling relation, e. g. massless QCD |72| or a fundamentally charged 
massless scalar coupled to Yang-Mills theory in Landau gauge |85l I86j . but this does not apply for the MAG. In the 
cited examples this is only possible because the matter fields do not couple to the ghost directly. For i5c = the same 
argument applies, which leaves as only possibility Sa + Sb ^ and Sa + Sc — leading to Saabb = Saacc — 0. Here 
the fields in the subscript determine the legs of the vertex to which the IRE belongs. It remains to fix the IREs of the 
other vertices. For vertices with an even number of legs a unique solution exists that corresponds to the maximally 
IR divergent solution. How this is realized is shown in fig.js] For a diagram with n^/ns/nc (all even) legs of type 
A/B/c one adds the necessary number of A ABB and AAcc vertices to a sunset diagram with a bare AABB or AAcc 
vertex. Since the added vertices do not scale, the additional internal propagators give exactly the necessary IREs to 
match the maximally IR divergent solution. 




FIG. 5: Consecutively inserting pairs of A, B and c fields shows that for graphs with an even number of legs the maximally IR 
divergent solution is realized. 

The situation is more intricate for vertices with an odd number of legs, starting already with ABB and Acc, which 
have two possible solutions: The one with the IRE — k/2, where k := 6b = Sc, corresponds to the maximally IR 
divergent solution. The other one has an IRE for the vertex functions ABB and Acc. This duality continues for 
higher vertex functions, where one can always show that there is an upper bound on the IRE equal to the lower bound 
plus k/2. The reason is purely combinatoric: The definitely leading graphs that can be identified unambiguously 
contain the AABB and AAcc vertices, but it is not possible to construct out of these a diagram with an odd number 
of legs. Thus, at the end there has to be always at least one vertex with an odd number of legs, which makes the 
analysis ambiguous. 

The IR scaling solution for SU{2) can be summarized as follows: 

• —Sa = Sb = Sc ■— k > 0, 

• d{nA, riB, ric) = \{nA - ub - nc)n {ua even), 

• 6{nA, riB, ric) = \{nA -ub - nc + ri)K (ua odd), 

where rj can be either or 1 and d{nA,nB,nc) is the IRE of a vertex function with ha diagonal gluon legs, ub 
off-diagonal gluon legs and ric ghost legs. We introduced the quantity k > in analogy to the usual notation in 
Landau gauge. For SU(2) the B and c fields only appear in even numbers. Since the IR leading behavior of DSEs is 
determined by four-point vertices, the dominant contributions arise from two-loop diagrams, namely the sunset and 
the squint diagrams with a bare AABB or AAcc vertex. However, there is the possibility that the squint diagrams are 
IR subleading. This can be traced back to the ambiguity of the vertices with an odd number of legs. Nevertheless the 
scaling relations Sa + 6b — and Sa + Sc = are the same and determined by the sunset diagrams alone. The question 
is: Are the squints IR leading as well or IR subleading due to the appearance of two three-point functions? The latter 
possibility is an example for a solution where the lower bounds for the IREs are not saturated. It is interesting that 
in RGEs the leading diagrams in propagator equations are tadpole diagrams in this case. 

To understand the ambiguity of vertices with an odd number of legs better we have a closer look at the DSE for the 
ABB vertex. It turns out that there are only three types of possibly leading diagrams: Those containing five-point 
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functions (directly coming from the leading sunsets of the propagator equations), those being reduced to Sabb upon 
inserting all results known up to now, and those being reduced to 56 abb + k. While the former two types are always 
leading, the latter diagrams only contribute to leading order if the maximally IR divergent solution is realized. On 
the other hand there are additional diagrams contributing to the leading order if Sab b = . The inequalities derived 
from these diagrams directly correspond to constraints already obtained above (as expected) and do not give any new 
information. The proof that the diagrams containing five-point functions are definitely leading is done via the DSE of 
the five-point function. As an example we explain it for the three-point diagram given in fig. [6] From the three-gluon 
vertex DSE we have 

Sabb < Sabbbb + k (40) 

and from the ABBBB vertex DSE 

5abbbb < Sabb — k Sabb > Sabbbb + (41) 

from which we can conclude Sabb = Sabbbb + The duality in the solution for vertices with an odd number of 
legs does not spoil the equivalence between Sabb and Sacc, since from the graphs in the corresponding DSEs given 
in fig.|7]one has both Sabb < Sacc and Sacc < Sabb, so that Sabb = Sacc- 




FIG. 6: One diagram of the DSE for the three-gluon vertex (left). Part of the DSE for the ABBBB vertex (right). 




FIG. 7: Diagrams of the ABB and Acc DSEs that allow to conclude Sabb ~ Sa, 



C. Inclusion of SU{N) interactions 

A special feature of the MAG is its gauge group dependence. It distinguishes it from Landau gauge, where there is 
no qualitative difference between SU{2) and SU{3) |62ffBl] . For the latter exists an additional structure constant in 
the MAG, which allows additional interaction terms in the Lagrangian. It has three off-diagonal indices, which clearly 
cannot be provided in SU{2), since there are only two off-diagonal fields. The step from SU{2) to SU{3) may thus 
change the results from the previous subsection. However, SU{3) already constitutes the most general case, since no 
other new interactions arise for higher SU (N) . 
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So take into account now the full plethora of interactions of SU(N) Yang-Mills theory. This leads to the following 
additional constraining inequalities for the propagator IREs: 

^Sb >0, ^b + 5, >0, 

\5a + ^Sb >0, + \5b + 5c >0. (42) 

They only contain combinations of fields that already occurred, e. g. + f can be written as combination of 
■^5a + Sb and Sb, of which we already know that they are non-negative. Therefore these inequalities are trivially 
fulfilled by the solution given above. Furthermore they cannot lead to any new solution, because a scaling solution 
from any of these terms, e. g. 1/2 (5^ -|- 1/2 (5^ + (^c = 0, would directly contradict the SU(2) inequalities. In the given 
example this would be 

Sa + Sc = {-Sb - 2Sc) + Sc = Sb - 5, > 0, (43) 

what directly contradicts the non-negativity of Sb and Sc, except both are 0. Thus for general SU{N) the same 
diagrams are IR leading as for SU{2). 

We established by now that for Green functions with an even number of legs only two-loop diagrams give the IR 
dominant contributions. This not only demands new methods for a numerical treatment, but also raises the number 
of dominant diagrams in higher vertex functions considerably, because the vast majority of diagrams possesses a bare 
four-point function. Already for the ghost-gluon vertex in SU (2) there are only four diagrams out of 27 including the 
tree-level vertex that are not contributing to the leading order if the maximally IR divergent solution is realized.^ 

We want to remark that the above solution is also valid if we would have started with the Lagrangian without 
the quartic ghost interaction, so that we would not have (5c > as a starting point. In this case one can analyze 
gluons only and establish that the two-loop diagrams with a bare AABB vertex are leading and Sa is non-positive. 
Including then the ghosts leads to the result that the AAcc vertex is leading and Sc is non-negative. Thus the ghost 
self-interaction is not decisive for the solution but makes the calculation easier, because right from the start it is 
known that Sc > 0. Nevertheless this remark is important, because the MAG is only realized by setting the gauge 
parameter a to and the bare quartic ghost interaction term is proportional to a. Thus some diagrams would vanish 
from the DSEs^, and to establish the inequality Sc> one needs the more involved analysis. 



D. Inclusion of the longitudinal dressing function of the B field 



In the analysis above we considered only one dressing function for the off-diagonal gluon propagator. Yet, if the 
gauge fixing parameter a is not set to zero (which is not possible in general due to terms proportional to 1/a in 
the gluonic vertices), the off-diagonal propagator has a longitudinal part. To take into account its own IRE Sbi we 
could add a fourth field and split the vertices accordingly into longitudinal and transverse parts. However, we went 
even one step further and inserted explicitly the expressions for the tree-level vertices, keeping the dressed vertices 
general. Projecting longitudinally and transversely, the number of transversal and longitudinal dressing functions is 
the same for each diagram, i. e. the equations for the two IREs Sb and Sbi are equal. This is in contrast to Landau 
gauge, where the bare three-gluon vertex is zero when contracted with three longitudinal projectors. In the MAG the 
three-gluon vertex has additional terms that do not vanish. Therefore we conclude that for the analyzed system the 
two dressing functions have the same IRE.^ 



EXISTENCE OF SCALING SOLUTIONS IN OTHER GAUGES 



The method developed in sec. |III| also allows to investigate other gauges. Again we disregard the possibility of 
potential cancellations and neglect bare two-point functions in the DSEs. 



® This applies only for that version of the ghost-gluon DSE, where the first differentiation is done with respect to the diagonal gluon. 
Note that terms proportional to 1/a exist, but they do not appear in the relevant cases to make contributions including a bare quartic 
ghost vertex. 

* This point could be invalidated in case contributions vanish nevertheless when the full (unknown) vertices are inserted into the DSEs. 
This would imply some yet unknown symmetry of vertex functions not taken into account here. 
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p ^= ^ ^ p 



FIG. 8: The ghost-gluon vertex DSE in Landau gauge. 



The Landau gauge was already used above as an example to show how the method developed in sec. |III| works and 
its unique scaling solution |18| [19] was established. We would like to add here an observation concerning the IRE 5gg 



of the ghost-gluon vertex: Using constraints from eq. (15 



25gg + 26gh + 6 gl > 0, 5gg + Ssg + 6 gh + 25 gl > 0, (44) 

in the ghost-gluon vertex DSE where the bare vertices are attached to an external ghost leg, see fig. |8] one can show 
that the contributions from the two triangle graphs yield an IRE greater or equal to zero. A similar constraint derived 



from eq. (15) states that the ghost-gluon scattering kernel cannot be so IR enhanced as to render the ghost-gluon 



vertex IR divergent: 

^2gh2g + Sgl + 5gh > 0. (45) 

Thus there is no contribution that is more IR divergent than the bare vertex and the ghost-gluon vertex dressing is 
IR finite. 

A natural extension is the step to general linear covariant gauges away from Landau gauge. In this case the gluon 
propagator has an additional longitudinal part. Due to gauge invariance this part acquires no dressing and stays 
bare [Tj. However, one can implement the longitudinal part into the presented formalism by a new IRE Sgij. Since 
the bare ghost-gluon vertex has a non-vanishing longitudinal part one gets an inequality connecting the longitudinal 
gluon IRE and the ghost IRE: 

ISgu + Sgh > 0. (46) 

The absence of a dressing function for the longitudinal part can be taken into account by setting Sgi^i to zero. Thus 
one has 

Sgh > (47) 



and no scaling relation can be found except the trivial one Sgh = 6g = 0. If inequ. (46 1 is indeed valid, the only way 
to circumvent this is to allow the longitudinal part to get dressed. This leads to Ogi— 6gij ~ —^Sgh, a result that 
was already found in ref. [87^ . 

As a final example we consider ghost anti-ghost symmetric gauges |88| [55] . In this gauge the longitudinal part of the 
gluon propagator gets dressed, but there is an additional quartic ghost interaction that again prohibits the existence 
of a scaling relation, because it leads to the non-negativity of the ghost IRE. As the two gluon IREs are non-negative, 
too, the only consistent solution is Sgi — 6gi^i — 6gh — 0. This result was already reported in ref. [S7], where bare 
four-point vertices were employed for the analysis. Our argument shows that without cancellations a different dressing 
of the vertices would not cure the failing of a scaling solution. 

These two examples allow two possible conclusions: Either the existence of a scaling solution is a special feature of 
some gauges, or there exist indeed reasons why the constraining inequalities do not hold, e. g. due to cancellations or 
the infiuence of kinematic divergences. In this context it seems promising to pursue interpolating gauges, which allow 
to investigate how confinement mechanisms in different gauges are linked. Speaking in terms of IREs the interactions 
that are switched on and off give the relevant information. However, one cannot expect from this point of view that 
there is a continuous transition, since a non- vanishing interaction always leads to the same scaling relations regardless 
of the numerical size of its coefficient, which is changed by the parameters interpolating between different gauges. An 



example is interpolation between Landau gauge and the MAG, see eq. (36) or ref. [55 . The leading vertices in the 
MAG are the AABB and AAcc vertices. While the two-ghost-two-gluon interaction is absent in Landau gauge at 
tree-level, the combination of all four-gluon interactions {AABB and BBBB) transforms to the standard four-gluon 
vertex of Landau gauge. Therefore the solution is expected to jump from IR enhanced diagonal gluons to IR enhanced 



ghosts as one sets the interpolation parameter Q in eq. (36) to zero. Whether the picture is really so simple will be 
investigated in the future. 
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VI. CONCLUSIONS 

We presented the IR scaling solution for the MAG, which gives a unique qualitative result for the propagators and 
vertices with an even number of legs. Vertices with an odd number allow several solutions which do not influence 
the general picture of the IR region, namely that the diagonal gluon is the dominant degree of freedom similar to the 
ghost propagator in Landau gauge. This behavior is in agreement with the hypothesis of Abelian dominance |51J, 
which states that the Abelian degrees of freedom should dominate in the IR regime. A necessary condition for the 
existence of a scaling solution is the cancellation of the bare propagator in the diagonal gluon propagator DSE. We 
speculated that this might be realized by a similar mechanism as the horizon condition in Landau gauge. 

Our analysis is to our knowledge the first investigation of the IR regime of the MAG for the physical gauge group 
SU{3). However, we did not find a qualitative difference to the simpler case of SU{2), which was not a priori clear 
due to the higher number of interaction terms in the Lagrangian for SU{3). These terms do not change the scaling 
solution obtained for SU{2) but only extend it to the new vertices. Furthermore they do not allow for additional 
solutions. So the qualitative results are obtained for general SU (N) and the emerging structure of DSEs in the IR is 
the following: The leading diagrams have two loops and a bare vertex that connects two diagonal and two off-diagonal 
fields. This implies that in general the major part of the diagrams in a DSE contributes to the leading IR order, since 
with the number of external legs the number of two-loop diagrams grows faster than that of one-loop diagrams. 

The method we used is based on the IREs of propagators and vertices and reduces an infinitely large set of DSEs to 
the relevant information in form of inequalities for the IREs. Thereby we proved that there is a lower bound for the 
IREs of vertices determined solely by the type and number of legs of the vertices. We called this bound the maximally 
IR divergent solution. We also showed that a scaling solution always leads to scaling relations linearly connecting the 
IREs of the propagators, which tells us that there always has to be at least one IR enhanced propagator. However, 
under certain circumstances (not applying to the MAG) it is possible that only some propagator IREs fulfill a scaling 
relation and part of the theory remains trivial in the IR. We want to stress that in the whole analysis we did not 
truncate the system in any way, but made only two assumptions: The one about the cancellation of the bare diagonal 
gluon propagator and that there are no cancellations for IR leading terms. 

It is also worth mentioning that the method suggested in the present article has been extended successfully to the 
case of Lagrangians which have a mixing of fields at the level of two-point functions [T5] as it happens, for example, 
in the Gribov-Zwanziger action [90 . This action realizes an improved gauge fixing, because the integration in field 
configuration space is effectively limited to the first Gribov region. It could be shown that this restriction does not 
interfere with the Dyson-Schwinger analysis [H], as already suggested earlier by Zwanziger [H], and the obtained 
solution agrees with the usual result from Faddeev-Popov theory [IHSl [TT] . 

An important point is that our analysis is only concerned with the IR part of the theory. In order to prove that 
the found solution really exists, it has to be shown in a numerical treatment over the whole momentum regime that it 
connects to the intermediate and ultraviolet momentum regime. But it is already a remarkable result that we could 
find a non-trivial IR scaling solution. Of course there may be other solutions besides an IR scaling solution as well, 
like for example the decoupling type of solution found by the refined Gribov-Zwanziger framework |58| and lattice 
simulations [54]. This issue might be connected with the choice of boundary conditions for the functional equations. 
In the MAG the boundary condition could be implemented most likely in the DSE of the diagonal gluon propagator 
by choosing a value for the dressing function at zero momentum. Also the possible presence of additional singularities 
in different kinematics has not been investigated. 

The next step would be of course the numerical solution of the equations. This presents due to the involved 
dynamics in the MAG a considerable challenge and is therefore left for future studies, which will have as a basis the 
results obtained here. The fact that two-loop diagrams are important in the IR prohibits the straightforward use of 
known non-perturbative renormalization methods. It constitutes an unexpected obstacle that is equally or even more 
complex than the high number of terms in the equations. 

We also applied the presented method to other gauges, where it cannot be a priori clear whether a scaling solution 
exists. Indeed our examples (linear covariant gauges, ghost anti- ghost symmetric gauges) are such, that no scaling 
solution can be found without further improvements to the employed technique. If thus a scaling solution is specific 
only to some gauges or if a more detailed analysis would reveal some yet unknown mechanism to allow for a scaling 
solution also in these gauges has still to be determined. Interpolating gauges might be an interesting application in 
this context. 
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Appendix A: Derivation of the IRE formula for an arbitrary diagram 



We derive the formula for the IRE Sy of a generic diagram v given in eq. ( 13 1 . We start by adding up the exponents 
of all momenta contributing in an /-loop diagram, i. e. the integral measures give a contribution ld/2, dressed 
propagators {Si — 1) and dressed vertices Si-^,,,i^', canonical dimensions are taken into account by Ci-^,,,i^ — 2 — 

i vertices, r'>3 

+ '^il...lr^il...t. ~ Cy. (Al) 

vertices ,r>3 

Hi are the number of internal propagators with IREs Si, whereas the numbers of vertices are 71*1 Superscripts d 
and b stand for dressed and bare, respectively. In case none is given, we refer to both. The sums '^^^^ti^f.^ r>3 extend 
over all vertices with at most r legs; mi being the number of external legs, the maximal value is + 2, since the 

highest vertex functions in DSEs and RGEs have at most two legs more than the vertex function on the left-hand 
side of the equation. 

We will use the following topological relations to rewrite the equation above: The number of vertices and propagators 
of a certain field (j)i are connected by 

li E , (A2) 



2 



1 vertices ^r'>?> 



where - denotes the number of times the field (pi appears in the vertex . . . c^^^. The number of loops can be 
replaced using 

l^'^Ui + l- E (A3) 



vertices ^r~>'3> 



vertices,r>3 



d 



Plugging these expressions into eq. (Al) yields 

^"^IEM E k^-'^n,,...,^-mi\+l- E 

\ 2 \ vertices, r'y^ I 

+ E M E fcr"''"^!..-. - ^^ - 1)+ 

i \ vertices ,r'> 3 j 

+ E <...v(v.v + 2-0+ E <...v(2-D-2 + ^E-^ 

vertices ^r> 3 vertices, r'>3 i 

- 2) (1 - ^ E - ^ E + E <..... (-1 + - + 2 - i 

\ i / i vertices,r>3 



+eU E + E <...J-^+2-0. (A4) 

i \vertices,r>3 J vertices^r>3 

Using J2i kl^'"^'' = r leads to 



Sr = {i-2 
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E 

vertices,r'>3 



"n.-.v I 1 7-1) (^-2) + ,5,,...,^ + 



E <...J(i-i)(-2) + iE^r^^4 

vertices.r>3 \ i / 



(A5) 



In four dimensions this result reduces to eq. ( 13 I 



Appendix B: Constraints from RGEs 

We outline the derivation of the constraints for four- and higher n-point functions and include the dimension d 
explicitly. The resulting expressions match exactly the coefficients of nf_^ in eq. (13). 
For four-point functions we get from fig. [9] the following two inequalities: 

i~'2]<6AABB. (Bl) 



Combining them yields 



2 

The inequalities for the three- and four-point functions fulfill 



5aabb < ^Sabcd + (5c + (^zj + ^ - 2. (B2) 



Sabcd + 1{5a + 5b + 5c + Sd) + 2(^~1]> 0. (B3) 



1) ('^ - 2) + <5,,...,„ + ^'"'^^^ ^ 0, (B4) 



d 

4 V ' ' 2 

3 

where r = 3, 4 denotes the number of legs the n-point function has. Higher n-point function also obey this inequality 
as can be shown by induction. 

For the proof we need the following two inequalities which can be inferred from fig. |10| 

Sabcd- - < Sabcdee--- + + - - 2, (B5) 
Saabbcc--- < 26abcde -- + 6d + Se + ^ - '2. (B6) 



The dots represent further legs as indicated in fig. 10 Note that these inequalities are generalizations of eqs. (Bll 
and (B2|. The first inequality can be used to write down an equation for the vertex AABBCC . . . as appearing in 
the second inequality: 

Saabb--- — Sc — — "^j < Saabbcc--- < '2Sabcde--- + Sd + Se + ^ — "2, (B7) 




FIG. 9: Parts of the FRGEs of generic four-point functions. 
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n 



D 




+ ... 

D 

- C b) 
FIG. 10: Parts of the FRGEs of generic n-point functions. 





AABB- 



2 < Sabcde - 



(B8) 



This inequality connects an n-point function with a (2n-6)-point function. The goal is to rewrite the equation such 



that only propagator IREs and the IRE of the n-point function remain. For this one can successively use eq. (B5 1 to 
replace the IRE of the remaining other vertex and arrives at 



1 



{Sa + Sb + Sc + Sd + Se + ■ ■ ■) - i -r - l] (r - 2) < Sabcde 



(B9) 



where r gives the number of legs of the n-point function. The coefficient of d/4 — 1 is obtained as follows: One 
gets 2((i/4 — 1) in eq. (B8) and d/4 — 1 from every further application of eq. ( |B5[ ). The latter has to be done until 
2ri — 6 — 2a; = 2, where x is the number of iterations. Hence the total coefficient of o?/4 — 1 is 2 + x — n — 2. This 



establishes eq. (B4| for all n-point functions. 



Appendix C: Inequality from the leading diagram 



We give here the derivation of eq. (pl| . For the leading diagram in a propagator DSE we get from eq. (|A5|) 



vertices \ i / dressed 



(Cl) 



Using rrii — 2 and miSi — 2Si this can be written as 



vertices \ i / dressed 



(C2) 



This is a general relation for the leading diagram, whichever this might be. Now we express the IREs of the vertices, 
by their lower bound, derived from eq. (A5l, to arrive at a new inequality: 



vertices \ i / 

E <....f-^EW-^^- 

dressed \ j 

vertices 



(C3) 



The right-hand side depends on dressed vertices only, indicated by the d superscript of n. On the other hand the 
left-hand side sums over dressed and bare vertices, so that in total only the bare vertex remains: 



> 0, 



(C4) 
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where we used J^i k^'"'^'' = r. For d = A this is exactly eq. (21 1 
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